Boundary-induced nonequilibrium phase transition into an absorbing state 
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We demonstrate that absorbing phase transitions in one dimension may be induced by the dy- 
namics of a single site. As an example we consider a one-dimensional model of diffusing particles, 
where a single site at the boundary evolves according to the dynamics of a contact process. As 
the rate for offspring production at this site is varied, the model exhibits a phase transition from 
a fluctuating active phase into an absorbing state. The universal properties of the transition are 
analyzed by numerical simulations and approximation techniques. 



Nonequilibrium phase transitions differ significantly 
from ordinary transitions at thermal equilibrium. For 
instance, under non-equilibrium conditions continuous 
phase transitions may occur even in one-dimensional sys- 
tems. A well-known example is the contact process for 
epidemic spreading 1] , where diffusing particles multiply 
at rate A and self-annihilate at rate 1. Depending on A, 
the contact process is either able to sustain a positive sta- 
tionary density of particles or it approaches a so-called 
absorbing state without particles from where it cannot es- 
cape. The active and the absorbing phase are separated 
by a continuous transition belonging to the universality 
class of directed percolation (DP) [3, H, 0], which plays 
a paradigmatic role like the Ising model in equilibrium 
statistical mechanics. Recently, the critical behavior of 
DP was confirmed experimentally for the first time by 
Takeuchi et al. [Bj]. 

As continuous phase transitions involve long-range cor- 
relations, boundary effects may play an important role. 
In the context of absorbing phase transitions previous 
studies focused primarily on DP confined to parabo- 
las 0, J?J , active walls Q , as well as absorbing walls and 
edges [9l [lO( . Although such boundaries influence the dy- 
namics deep in the bulk, the universality class of the bulk 
transition is not changed inherently, rather it is extended 
by an additional exponent describing the order parame- 
ter near the boundary. A completely different situation 
is encountered in systems where boundary effects induce 
a new type transition which would be absent without the 
boundary [111 ] . Such boundary-induced phase transitions 
have been studied for example in models for diffusive 
transport [13, El and traffic flow Q. 

In this Letter we present an example of a boundary- 
induced phase transition from a fluctuating phase into 
an absorbing state. To this end we consider a simple 
one-dimensional model, where the leftmost site evolves 
in the same way as in the contact process while particles 
in the bulk diffuse according to a symmetric exclusion 
process. Varying the rate for offspring production at the 
leftmost site the model exhibits a non-equilibrium phase 
transition from a fluctuating active phase into an absorb- 
ing state with a non-trivial critical behavior. A similar 
problem with catalytic creation and pair-annihilation at 
a single site in the center and diffusion in the bulk was 



studied in [l5| by field-theoretic methods. 

Definition of the model: The model is defined on an 
semi-infinite one-dimensional chain of sites i = 0, 1, 2, . . . 
which are either empty (s, = 0) or occupied by a particle 
(si = 1) (see Fig. Q] ). Starting at time t = with a 
single particle at the origin (si = J^o) the model evolves 
by random-sequential updates as follows. For each up- 
date one of the particles is randomly selected. If the 
selected particle is located at the leftmost site i = 0, it 
undergoes the same dynamics as in a standard contact 
process, namely: 

(a) With probability p = A/(l + A) a new particle is 
created at the right neighbor, provided that this 
site is empty. This can be done by setting Si := 1. 

(b) Otherwise, the particle at the leftmost site is de- 
stroyed by setting sq := 0. 

Else, if the selected particle is not located at the origin, it 
diffuses according to a symmetric exclusion process, i.e., 
it jumps to a randomly chosen nearest neighbor, pro- 
vided that the target site is empty. As usual in models 
with random-sequential dynamics, each attempted up- 
date corresponds to a time increment of 1/N(t), where 
N(t) is the actual number of particles. On a computer 
the dynamical rule defined above can be implemented ef- 
ficiently by using a dynamically generated list of particle 
coordinates, eliminating possible finite-size effects. 

Phenomenological properties: In the bulk the sym- 
metric exclusion process preserves the number of par- 
ticles, whereas this conservation law is violated at the 
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FIG. 1: Definition of the model: Particles diffuse on a semi- 
infinite one-dimensional chain according to a symmetric ex- 
clusion process. The only exception is the leftmost site, where 
particles multiply and annihilate as in a contact process. 
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FIG. 2: Typical temporal evolution of the model at criticality, 
starting with a single particle at the leftmost site. The color 
scale visualizes the particle density. At the leftmost site one 
observes intermittent bursts of activity. 



boundary, where offspring production and removal com- 
pete one another. For p = the leftmost site acts as 
a sink where particles disappear, thereby depleting the 
whole system diffusively until the dynamics reaches the 
absorbing state without particles. On the other hand, for 
p = 1 the leftmost site is permanently occupied, provid- 
ing a steady source of particles at the left boundary so 
that the system approaches a fully occupied stationary 
state. In between it turns out that the (infinite) system 
is able to maintain a non-vanishing stationary density of 
particles even for finite values of p down to a well-defined 
critical threshold p c . 

Starting with a single particle at the boundary, the 
process evolves as follows. Initially the particle at the 
leftmost site either disappears or it creates another par- 
ticle at its right neighbor. As soon as this freshly created 
particle diffuses away into the bulk, the particle at the 
leftmost site may create and send out further particles 
until it disappears by spontaneous removal. The average 
number of newly created particles depends on p and is of 
order 1 in the stationary state. 

Each created particle performs a one-dimensional ran- 
dom walk in the bulk, which in one dimension is bound 
to return to the origin after finite time. The returning 
particles may either disappear or release another bunch 
of particles. As demonstrated in Fig. [21 particles are not 
created continuously but in form of intermittent bursts. 
Apparently these irregular bursts are responsible for the 
nontrivial properties of the model. 

Seed simulations at criticality: The simplest order pa- 
rameter describing the phase transition is the occupation 
probability po(t) = (so(t)) of the leftmost site averaged 



over many independent runs. For small p« 1 this quan- 
tity is dominated by the first-return probability of a one- 
dimensional random walk which is known to decay with 
time as i -3 / 2 (see e.g. This power- law decay char- 

acterizes the inactive phase of the system. Contrarily, 
for large values of p, the returning particle is likely to 
multiply frequently, flooding the bulk of the system with 
freshly created particles and thereby maintaining a con- 
stant non-zero density. In between we find a phase tran- 
sition located at the critical point (see Fig. [3|) 

Pc = 0.74435(15), (1) 

at which po(t) decays as 

p (t)~t- a , a = 0.50(1), (2) 

suggesting the exact value a = 1/2. 

Another well-known order parameter is the survival 
probability P s (t) to find at least one particle in the entire 
system at time t. At the transition this quantity is found 
to decay algebraically as 

P s (t) = r 5 , 6 » 0.15(2). (3) 

This estimate shows a slight systematic drift and may be 
compatible with the rational value 6 = 1/6. 

Finally, it is useful to study the number of particles 
(N(t)) averaged over all runs. For small p, this quantity 
decreases as (N(t)) ~ t^ 1 / 2 , while for large values of p 
one finds an algebraic increase (N(t)) ~ t +1/2 because 
of the diffusive bulk dynamics. At the transition (N(t)) 
stays almost constant close to 1, suggesting a vanishing 
exponent 

(N(t)) ~ i e , 9 = 0. (4) 

Note that (N(i)) is averaged over all runs. If the number 
of particles was averaged over surviving runs it would 
actually grows as t 5 . 

We also determined the density profile p(x, t) in the 
bulk (see Fig. 0]). At criticality this profile turns out 
to be an almost perfect Gaussian distribution, obeying 
the scaling form p(x,t) = t~ x l 2 f(x 2 /t). This indicates 
a simple diffusive behavior in the bulk. On the other 
hand, the density of neighboring pairs of particles is not 
a Gaussian, demonstrating that the random walks are 
mutually correlated. 

Off-critical seed simulations: As can be seen in Fig. [3l 
for p > p c the average density at the boundary first de- 
creases algebraically, goes through a minimum, then in- 
creases again until it reaches a stationary value. Surpris- 
ingly, the time at which the minimum is reached scales 
roughly as (p— p c ) ~ 3 while the stationary value is reached 
at a typical time that scales as (p — p c )~ A ■ Therefore, 
it is impossible to produce a data collapse by plotting 
p(t) I (p — pc) 13 versus t{p — p c ) u ^ . However, collapsing the 
crossover from increase to saturation, one would consis- 
tently get the exponents (3 = 1 and = 4. 



3 




t t t 

FIG. 3: Numerical simulation of the process starting with a single active site at the origin for different values of A = p—p c - The 
graphs show the density at the leftmost site po(t), the survival probability P„(t), and the number of particles (N(t)) averaged 
over all runs. 




FIG. 4: Left: Data collapse of the rescaled profiles of the 
particle density at criticality for to = 64, 128, . . . , 8192 (blue) 
compared to a Gaussian distribution (red). Inset: The same 
data collapse in a double-logarithmic representation. Right: 
Density of particles (blue) and pairs (green) at to = 10 6 , 
demonstrating the presence of correlations which decay in 
space as x _1,/2 , indicating that /3/f± = 1/2. 



Homogeneous initial state: Starting with a fully oc- 
cupied lattice at criticality, the special dynamics at the 
leftmost site gradually depletes the system, leading to a 
slow decay of the particle density at the boundary. In 
numerical simulations we find that the density decays 
slowly as t~ s , where S = 0.15(2) is the same exponent 
as in Eq. ([3]) which describes the survival of a cluster 
generated from a single seed. As will be explained in 
a forthcoming publication, this can be traced back to a 
duality of the two situations under time reversal. 

Mean field analysis: In a simple mean field approx- 
imation the n-site probability distribution is approxi- 
mated by the product of n single-site probabilities, ne- 
glecting correlations. Defining r/i(t) as the first moment 
of the probability distribution at site i, the mean field 
equations read 

^ = -(l-p)*?o + ^i(l-%), (5) 
= P»7o(l-»7i) + ^(»72+»7o»7i-2?7i), (6) 
= ^ (Vi+i + Vi-i ~ 2j 7<) for * = 2, 3, (7) 



Solving these equations for 77^ (0) = 8i t o we find the crit- 
ical point is p c MF = 1/2, where equation ([6]) reduces to 
a diffusion equation, reproducing the critical exponents 
_ an( j q mf _ \J2. However, starting with a fully 
occupied lattice 77^ (0) = 1 one gets a decay r]o(t) ~ i~ 1//4 , 
differing from the simulation result [l8j |. 

Possible relation to a non-Markovian process: To un- 
derstand the transition of the model from a different 
point of view, let us now adopt the perspective of the 
leftmost site. If this site is active, it may create new 
particles, sending them out for random walk in the bulk. 
From the prospect of the leftmost site the specific tra- 
jectory of this random walk does not matter, the only 
question of interest will be at which time the particle 
returns to the origin. 

Let us now assume that the diffusing particles in the 
bulk do not interact. For a symmetric exclusion process 
this approximation is justified if the particle densities are 
sufficiently small. With this approximation a particle 
emitted at the leftmost site will return after a time At 
which is distributed algebraically as 

P(At) ~ (At)" 3 / 2 . (8) 

Following [l5| the problem can be reformulated as a 
single-site process with a non-Markovian dynamics. Let 
s(t) = 0, 1 denote the occupancy of a single site at time 
t £ N, which can be implemented as a one-dimensional 
array s [t] on a computer. The array is initialized by 
s(t) :— St t o, corresponding to a single particle at the 
boundary. The single-site model then evolves according 
to the following dynamical rules: 

1. Select the lowest t for which s(t) = 1. 

2. With probability /i generate a waiting time At ac- 
cording to the distribution ([8]), truncate it to an 
integer, and set s(t + At) := 1. 

3. Otherwise (with probability 1 — fi) set s(t) := 0. 

These steps are repeated until the system enters the ab- 
sorbing state or t exceeds a predetermined maximal time. 
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Simulating this non-Markovian single-site process us- 
ing a dynamically generated list we can go up to 10 12 
time steps, finding the critical point \i c — 0.57426(1), 
the correct exponent a = 0.500(5), as well as a consis- 
tent exponent for the survival probability 5 = 0.16(1). 
Moreover, the off-critical properties of the original model 
are faithfully reproduced. This suggests that the non- 
Markovian process defined above may be even equivalent 
to the original model regarding its asymptotic critical 
behavior. This is surprising since the approximation ig- 
nores the exclusion principle of the random walkers in 
the bulk. 

Relation to a non-Markovian Langevin equation: Let 
us finally describe the single-site process in the contin- 



uum limit. As shown in previous studies (see e.g. 17 1 
and references therein), a non-Markovian dynamics by 
algebraically distributed waiting times -P(At) ~ At _1_K 
is generated by so-called fractional derivatives <9 t K which 
are defined by 



dtP(t) = 



1 



Mk) 



dt'r 



■>[p(t)-p(t-t')}, (9) 



where n £ [0, 1] and A/ii(/t) = — T(— k) is a normalization 
constant. This suggests that the non-Markovian single- 
site model may be effectively described by a DP-likc 
Langevin equation without space dependence, in which 
the local time derivative is replaced by a fractional deriva- 
tive with At = 1/2 generating temporal Levy flights: 



di /2 P (t) = ap(t)-p(ty + at)' 



(10) 



Here the parameter a plays the role of /i — fi c , the second 
term accounts for the fact that the leftmost site cannot be 
activated twice, and £ is a multiplicative noise with cor- 
relations (£,(t)^(t')} — p(t)8(t — t'). Dimensional analysis 
confirms that the noise is relevant under temporal rescal- 
ing, supporting the expectation that the model exhibits 
a non-mean-field properties. We note that this Langevin 
equation can be converted by standard techniques into a 
Fokker Planck equation of the form 



dl' 2 P{p,t) 



(ap - p 2 )P(p,t) 



2 " 



pP(p,t) 



This equation can be decoupled by a Laplace transforma- 
tion in t but so far we were not able to solve the resulting 
ordinary differential equations analytically. 

Towards a scaling picture: Based on these results we 
conjecture that the full model is described by an order 
parameter field at the boundary with the exponent 0=1 
and a response field with the same exponent j3' = 1. The 
diffusive dynamics in the bulk, characterized by a dynam- 
ical exponent z = 2, is slaved to the boundary dynamics 
and thus it does not induce additional order parameter 
exponents. Moreover, the numerical results indicate that 
the scaling exponents are given by v\\ = 4 and v±_ = 2. 



The decay exponent a describes a two-point function and 
hence it is given by a — ((3 + [3')/v\\ = 1/2, while the slip 
exponent is consistently described by the hyperscaling 
relation 9 = ((3 + /3' — =0. The only exponent, 

which is so far not explained within this picture, is the 
survival exponent 8 « 0.16. We believe that this can be 
traced back to the non-Markovian character of Eq. (| 10[) , 
which requires a new understanding of the survival prob- 
ability. 

To summarize, we have studied a model that exhibits 
a novel class of boundary-induced phase transition from 
a fluctuating phase into an absorbing state. This is prob- 
ably the simplest non-trivial absorbing phase transition, 
much simpler than ordinary DP, but nevertheless exhibit- 
ing properties which cannot be explained within mean 
field theory. 
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